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OC I We first refine the fixed concept in the hterature that the usage of the New- 

tonian potential in studies of black hole accretion is invalid and the general rel- 
Kj^ , ativistic effect must be considered. As our main results, we then show that the 

O ! energy released by neutrino annihilation in neutrino-dominated accretion flows 

is sufficient for gamma-ray bursts when the contribution from the optically thick 






"^ I region of the flow is included, and that in the optically thick region advection 

\Q I does not necessarily dominate over neutrino cooling because the advection factor 
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is relevant to the geometrical depth rather than the optical depth of the flow. 
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1. Introduction 

The fireball shock model (see, e.g., Meszaros 2002 and Zhang & Meszaros 2004 for 
reviews) has been widely accepted to interpret the gamma-ray and afterglow emitting of 
gamma-ray bursts (GRBs). Despite the successes of this phenomenological model, the central 
engine of the relativistic fireball is not yet well understood. Most popular models for the 
energy source of GRBs are in common invoking a hyperaccreting black hole. Accretion 
models in this context were first considered by Narayan, Paczyhski, & Piran (1992), and 
have been recently discussed by Popham, Woosley, & Fryer (1999, hereafter PWF), Narayan, 
Piran, & Kumar (2001), Kohri & Mineshige (2002), Di Matteo, Piran, & Narayan (2002, 
hereafter DPN), and Kohri, Narayan, & Piran (2005). 

PWF introduced the concept of neutrino-dominated accretion flows (NDAFs) and showed 
that the energy released by neutrino annihilation was adequate for GRBs. Their calcula- 
tions, however, were based on the assumption that the flow is optically thin for neutrinos. 
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As pointed out by themselves, this assumption breaks down for the mass accretion rate 
M > 10 Mq s~^. They mentioned that their estimate of the neutrino annihilation luminos- 
ity ~ 2 x 10^^ ergs s~^(see their Table 3) for M = 10 Mq s~^ should be taken as an upper 
limit, and the actual luminosity could be as much as a factor of 5 lower, i.e., ~ 4 x 10^^ 
ergs s~^. The NDAF model was reinvestigated by DPN, in which a bridging formula was 
adopted for calculating neutrino radiation in both the optically thin and thick cases. They 
showed that for M > 0.1 Mq s~^ there exists an optically thick inner region in the flow; 
and argued that for M > 1 Mq s~^ neutrinos are sufficiently trapped and energy advection 
becomes the dominant cooling mechanism, resulting in the maximum luminosity of neutrino 
annihilation which is only ~ 10^° ergs s~^(see their Fig. 6). Thus they claimed that the 
NDAF model cannot account for GRBs. 

How to understand the inconsistent results of PWF and DPN? We note that PWF 
worked in the relativistic Kerr geometry, but with the a priori assumption that neutrinos 
are optically thin; whereas DPN calculated the optical depth for neutrinos, but went back 
into the Newtonian potential and omitted totally the neutrino radiation from the optically 
thick region. The purpose of this Letter is to try to update partly the NDAF model. It 
is surely correct that the general relativistic effect must be considered in studies of black 
hole accretion, then we wish to know how important the effect of the neutrino opacity is in 
determining the luminosity of an NDAF. 



2. Assumptions and equations 

For simplicity, a steady state axisymmetric black hole accretion flow is considered as 
in PWF and DPN. We adopt that the general relativistic effect of the central black hole is 
simulated by the well-known pseudo-Newtonian potential introduced by Paczyhski & Wiita 
(1980, hereafter PW potential), i.e., $ = -GMbh/I-R - -Rg), where Mbh is the black hole 
mass, R is the radius, and R^ = 2GMbh/c^ is the Schwarzschild radius. Other assumptions 
about the flow are usual in the literature: the angular velocity is approximately Keplerian, 
i.e., Q = Qk = (GMbk/R)^ /{R — Rg); the vertical scale height of the flow is H = Cs/I^k, 
where Cg = (P/p) ^ is the isothermal sound speed, with P and p being the pressure and mass 
density, respectively; and the kinematic viscosity coefficient is expressed as z/ = aCsH, where 
a is the constant viscosity parameter. 

The basic equations describing the flow consist of the continuity, azimuthal momentum, 
and energy equations plus the equation of state. The continuity equation is 

M = -AnpHRv , (1) 
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where v is the radial velocity. With the assumption Vt = Qk, the azimuthal momentum 
equation is reduced to an algebraic form: 

TT 

V = -acs—r^g , (2) 

where / = l—j/{QR'^) and g = — din fix/din i?, with the integration constant j representing 
the specific angular momentum (per unit mass) accreted by the black hole. The equation of 
state is written as 

-^ ^ -'gas ~r Jrad + -^ dcg + -'i^ ; (,<jj 

where Pgas, -Prad, -Pdeg, and Pi, are the gas pressure from nucleons, radiation pressure of 
photons, degeneracy pressure of electrons, and radiation pressure of neutrinos, respectively. 
The energy equation is written as 

^; vis 'Vadv ~r (Vphoto ~r ^u ■ \^) 

The viscous heating Qvis and the advective cooling Qadv (for a half disk above or below the 
equator) are expressed as 

(5) 
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where T is the temperature, s is the specific entropy, Xnuc is the mass fraction of free 
nucleons, u^, is the neutrino energy density, and ^ oc — dlns/dlni? is assumed to be equal 
to 1 as in DPN. The quantity Qphoto is the cooling of the photodisintegration process, and 
Qu is the cooling of the neutrino radiation. We adopt a bridging formula for calculating Q^,, 
which is valid in both the optically thin and thick cases. Detailed expressions for Pgas, Prad, 
Pdcg, Pu, Qphoto, -^nuc, u^, and the bridging formula for Q^, are given in DPN. 

Equations (1-4) contain four independent unknowns p, T, H, and v as functions of R, 
which can be numerically solved with given constant parameters Mbh, M, a, and j, then all 
the other quantities can be obtained. In the following calculations we fix Mbh = 3Mq and 
a = 0.1. 



3. Invalidity of the usage of the Newtonian potential 

Most previous calculations for NDAFs (e.g., Narayan et al. 2001; DPN; Kohri et al. 
2005) adopted the Newtonian potential and did not take the integration constant j into 
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consideration. Kohri & Mineshige (2002) also used the Newtonian potential but considered 
j. Only PWF worked in the relativistic Kerr geometry as we mentioned already. In this 
section we refine the fixed concept in the literature, i.e., the invalidity of the usage of the 
Newtonian potential. We concentrate on three solutions corresponding to PW potential with 
j = l.ScRg (just a little less than the Keplerian angular momentum at the last stable orbit, 
^xlsflg = 1.837ci?g), the Newtonian potential with j = 1.2cRg (/xls/jg = 1.225c-Rg, cf. Kohri 
& Mineshige 2002), and the Newtonian potential with j = (DPN; Kohri et al. 2005), 
respectively. 

The variation of the optical depth r for neutrinos with R is drawn in Figure 1(a), for 
which the dimensionless mass accretion rate is m = M/{Mq s~^) = 1. The figure shows that 
the values of r in the Newtonian potential (the dotted and dashed lines) are significantly 
larger than those in PW potential (the solid line). The accretion fiow in PW potential is 
completely optically thin (r < 2/3), whereas for the Newtonian potential with j = there 
exists a wide optically thick (r > 2/3) region of i? < 15.4_Rg. We believe that the results with 
PW potential are more convincible since this potential is known to be a better description 
for a nonrotating black hole than the Newtonian potential. Our arguement can be further 
confirmed by Figure 1(b), which shows the variation of rj^ with rh, where r]i, = Lj^/Mc^ 
is the efficiency of energy release by neutrino radiation (before annihilation). As seen in 
the figure, rju in the Newtonian potential is much larger than that in PW potential. For 
j = 0, the former can reach a maximum value of 0.206 at rh = 0.45, which is far beyond the 
maximum possible efficiency in the Schwarzschild geometry {rj = 0.057) and is unphysical. 
In fact, by integrating the viscous heating Qvis from 3i?g to the infinite outer boundary of the 
fiow, we can obtain the theoretical maximum rji, for the above three solutions: 1/4 for the 
Newtonian potential with j = [from Eq. (14) of DPN], 1/12 for the Newtonian potential 
with j = 1.225ci?g [from Eq. (32) of Kohri & Mineshige 2002], and 1/16 for PW potential 
with j = 1.837cRg. Obviously the result with PW potential is the closest to the reality 
(0.057), while the results with the Newtonian potential are unreasonable. 

We conclude for the moment that the usage of the Newtonian potential is invalid in 
calculations for NDAFs at least at the following two points: (1) it would overestimate sub- 
stantially the optical depth for neutrinos; (2) it would lead to an unphysical efficiency of 
energy release by neutrino radiation. 



4. Effect of the optical depth on the neutrino annihilation luminosity 

Our method for calculating neutrino annihilation is similar to many previous works 
[e.g., Ruffert et al. 1997; PWF; Rosswog, Ramirez- Ruiz, & Davies 2003). Figure 2 shows 
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the variations of Ly (the upper thin sohd hne) and Lj^p (the lower thick sohd hne) with fn, 
where Ly is the luminosity of neutrino radiation before annihilation, and L^p is the luminosity 
of neutrino annihilation (which is the most important from the observational point of view), 
both of them are calculated with PW potential. The circles and triangles represent the 
results of PWF for Ly (empty) and Lyy (filled), respectively. It is seen that our results agree 
very well with that of PWF for m < 1, because PW potential is a good approximation for 
the Schwarzschild geometry. For m > 1, our results are lower than that of PWF. This is 
because they assumed neutrinos to be optically thin; while we use the bridging formula for 
Qu^ and there exists an optically thick region for fn > 1.2. According to our calculations 
Lyy varies from 3.9 x 10^*^ ergs s^^to 3.6 x 10^^ ergs s~^for 1 < m < 10, which implies that, 
based on the energy consideration, NDAF can indeed work as the central engine for GRBs. 
In particular, our L^y (3.6 x 10^^ ergs s^^) for m = 10 is in good agreement with PWF's 
"actual luminosity" (~ 4 x 10^^ ergs s""*^, as mentioned in § 1). 

For comparison. Figure 2 also shows Lyy in other three cases: (1) using PW potential but 
omitting the contribution from the optically thick region (the dotted line, the r > 2/3 region 
appears for m > 1.2); (2) using the Newtonian potential and including the contribution from 
the optically thick region (the dot-dashed line); and (3) using the Newtonian potential but 
omitting that contribution (the dashed line, the r > 2/3 region appears for m > 0.052). 
As known from § 3, the results of cases (2) and (3) are unreal because the usage of the 
Newtonian potential overestimates unphysically both r and Ly. It is seen that the omitting 
of the contribution from the r > 2/3 region reduces substantially Lyp in case (3), i.e., even 
with the overestimated Ly caused by the Newtonian potential; as well as in case (1), i.e., 
even the general relativistic effect is considered. This is probably the reason why DPN 
obtained Lyy in their Newtonian calculations which is insufficient for GRBs. We think that 
it is unfair to ignore totally the neutrino radiation from the r > 2/3 region. As DPN also 
stated, the neutrino emission is partially suppressed as the inner regions of the flow are 
becoming opaque. The trapping of neutrinos is a process that is strengthening gradually 
with increasing r, the value of r reaching 2/3 does not mean that all neutrinos are suddenly 
trapped, and the use of the bridging formula for Qy is exactly to calculate the neutrino 
radiation from both the optically thin and thick regions. A similar bridging formula has 
been widely used for calculating the radiation of photons in both the optically thin and 
thick cases (e.g., Narayan & Yi 1995). 
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5. Energy advection 

DPN argued that energy advection would become the dominant coohng mechanism when 
the flow is optically thick for neutrinos. As seen from their Fig. 3, however, it is not the case. 
For example, for rh = 1 the flow is optically thick at i? = 10-Rg, but the advection factor 
/adv = Qadv/Qvis at this radius is only ~ 0.1. In our opinion, whether cooling is dominated 
by advection or by radiation is not determined by the optical depth. For accretion flows in 
X-ray binaries and AGNs, it is known that /adv is relevant to the geometrical depth rather 
than the optical depth of the flow (Abramowicz, Lasota, & Xu 1986): 

/adv OC i^r . (7) 

Such a relationship can be well checked by the four representative types of accretion mod- 
els: the optically thick standard thin disk (Shakura & Sunyaev 1973) and the optically thin 
SLE disk (Shapiro, Lightman, & Eardley 1976) are both geometrically thin and radiation- 
dominated, i.e., /adv ~ 0; and the optically thick slim disk (Abramowicz et al. 1988) and 
the optically thin advection-dominated accretion flow (Narayan & Yi 1994) are both geo- 
metrically thick and advection-dominated, i.e., /adv ~ 1- We argue that this relationship 
should also work in the NDAF model with the following modiflcation in accordance with 
PW potential: 

/adv OC /^ = (f )'/"'^~' ' (8) 

where fn is called by us the geometrical depth factor, and f^^g^^ comes from the expression 
ofQvis [Eq. (5)]. 

As shown in Figure 3, the variation of /adv with R (the solid line) agrees very well with 
that of fff (the dot-dashed hne), but differs signiflcantly from that of r (the dotted line), 
clearly indicating that the strength of energy advection is relevant to the geometrical depth 
rather than the optical depth. It is also seen that, although the flow is optically thick for 
R < 25.9-Rg, advection dominates over neutrino cooling (/adv > fu, where fu = Qv/Qvis 
is drawn by the dashed line) only in a smaller region R < 5.1-Rg. Once again, this result 
supports the view that it is important to consider the role of the optically thick, but neutrino 
radiation-dominated region, e.g., 5.1i?g < R < 25.9-Rg in the example of Figure 3. 



6. Discussion 

We have shown that the usage of the Newtonian potential along with the omitting 
of neutrino radiation from the optically thick region would lead to unreal luminosities for 
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NDAFs, and that when the general relativistic effect is considered and the contribution from 
the optically thick region is included, NDAFs can work as the central engine for GRBs. 

In addition to its mass, a black hole may have its spin as the other fundamental property. 
We consider here only the nonrotating black hole, for which PW potential can work. PWF 
has shown that a spining (Kerr) black hole will enhance the efficiency of neutrino radiation. 
This strengthens our conclusion here that NDAFs into black holes can be the central engine 
for GRBs. 

We have tried to update partly the NDAF model by considering both the effects of 
general relativity and neutrino opacity. There are certainly other factors which may influ- 
ence the neutrino luminosity of an NDAF and we do not consider here, such as the electron 
degeneracy. We adopt a simple treatment for the electron degeneracy pressure in agreement 
with PWF and DPN. Kohri & Mineshige (2002) pointed out that it is important to include 
the effect of electron degeneracy that suppresses the neutrino cooling at high density and 
high temperature. Most recently, Kohri et al. (2005) considered the effects of both electron 
degeneracy and neutrino optical depth, and calculated the neutrino cooling, the electron 
pressure, and other physical quantities even in the delicate regime where the electron de- 
generacy is moderate; while in previous works as well as ours here the calculations can be 
made accurate only in the two opposite limits of extremely degenerate electrons and fully 
nondegenerate electrons. We wish to see in future studies how the electron degeneracy would 
affect our results here. 

We thank Ye-Fei Yuan for beneficial discussion and the referee for helpful comments. 
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Fig. 1. — Solutions in PW potential with the accreted specific angular momentum j = l.ScRg 
(the solid line), in the Newtonian potential with j = 1.2cRg (the dashed line), and in the 
Newtonian potential with j = (the dotted line), (a) Neutrino optical depth r as a function 
of radius R for the dimensionless mass accretion rate m = 1; (b) Efficiency of neutrino 
radiation rji, as a function of m. 
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Fig. 2. — Neutrino luminosity before annihilation L^, with PW potential (the thin solid line), 
neutrino annihilation luminosity L^p with PW potential and including the r > 2/3 region 
(the thick solid line), L^p with PW potential but omitting the r > 2/3 region (the dotted 
line), Lijp with the Newtonian potential and including the r > 2/3 region (the dot-dashed 
line), and Li^p with the Newtonian potential but omitting the r > 2/3 region (the dashed 
line) as functions of m. The empty circles and triangle denote L^ of PWF, and the filled 
circles and triangle denote L,yp of PWF. 
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Fig. 3. — Advection factor /adv, geometrical depth factor fn, neutrino cooling factor /j,, 
and r as functions of R for m = 5. The filled circle denotes the r = 2/3 position, i.e., 
R = 25.9R,. 



